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Abstract 



Qh ' We establish one-dimensional spectral, or "normal modes" , stability of ZND deto- 

•^f nations in the small heat release limit and the related high overdrive limit with heat 

release and activation energy held fixed, verifying numerical observations of Erpenbeck 
in the 1960s. The key technical points are a strategic rescaling of parameters convert- 
ing the infinite overdrive limit to a finite, regular perturbation problem, and a careful 
high-frequency analysis depending uniformly on model parameters. The latter recovers 
the important result of high-frequency stability established by Erpenbeck by somewhat 
different techniques. Notably, the techniques used here yield quantitative estimates well 
^- ' suited for numerical stability investigation. 

CO 

f^> ■ 1 Introduction 

in 

In this note, we establish one-dimensional spectral stability in the small heat release and 
high overdrive limits of ZND detonations, confirming numerical observations of Erpenbeck 
going back to [Er2], but up to now not rigorously verified. In the process, we recover by a 
somewhat different argument the fundamental result of Erpenbeck [Er3] that detonations 
are one-dimensionally stable with respect to sufficiently high frequencies. 

The basic argument for stability in the high overdrive limit, based on a strategic 
rescaling of parameters converting the problem to a small-heat release/small activation 
energy/maximal shock strength limit on a bounded parameter range, was indicated already 
in [Zl]. However, the result seems sufficiently fundamental to warrant an exposition giving 
full detail. In particular, the discussion of [Zl] concerned only bounded frequencies, for 
which stability follows by a simple continuity argument. Stability for high frequencies can 
be concluded from a well-known result of [Er3], which, restricted to the one-dimensional 
setting, asserts that instabilities cannot occur outside a sufficiently large ball. However, 
the arguments of [Er3], based on semiclassical limit and turning point theory, 1 do not 
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1 Specifically, validation of a formal asymptotic expansion in one over the frequency for complexified x. 



readily yield quantitative estimates on rates of convergence or dependence on parameters. 
As continuous dependence on parameters of the radius outside which instabilities are ex- 
cluded is crucial for the limiting argument described above, it seems useful to revisit the 
high-frequency limit problem in greater detail. 

Moreover, as pointed out in [CJLW], the issue of matching at x — > — oo of the formal 
asymptotic solution with the solution prescribed by the required behavior at spatial infinity 
appears to require a bit more discussion beyond what is given in [Er3], where it is concluded 
simply from the observation that the limits as x or frequency go to infinity commute. This 
argument seems to require either uniform convergence on all of x £ (— oo,0] of the formal 
asymptotic series as frequency goes to infinity, which (to us) does not appear obvious, or 
else uniform estimates independent of frequency on behavior as x — > — oo. 

These aspects (uniform dependence and uniform convergence on (— oo,0]) of high- 
frequency behavior are the main issues addressed here, where they are treated by a careful 
application of the asymptotic ODE techniques developed in [GZ, MeZl, MaZ3, PZ, Zl]. 
These in turn are natural outgrowths of the classical asymptotic ODE techniques devel- 
oped by Coddington, Levinson, Coppel, and others, as described in [CL, Co] and references 
therein, including those cited by Erpenbeck [Erl, Er2] in describing behavior as x — > — oo. 2 
Our arguments are qualitatively different from the ones of Erpenbeck based on semiclassi- 
cal limit/turning point theory, making use at a key point of the exponential convergence 
of profiles to a limit as x — > — oo. Specifically, in the neutral case that diagonal elements 
have uniformly small spectral gap, we apply a variable-coefficient version (Lemma 4.3) of 
the conjugation lemma of [MeZl] to close the argument, extending and refining the related 
constant-coefficient gap lemma estimates of Proposition 5.7, [CJLW], applying to a single 
mode. 3 

Notably, the simple and concrete estimates thus derived yield quantitative bounds of 
potential use for numerical stability investigations. Recall [HuZl] that computational in- 
tensity of numerical stability computations increases rapidly with increasing frequency, so 
that bounds on frequency are of considerable practical interest for applications. See [BZ] 
for a first effort in this direction in the simplified context of Majda's model [M]. 

2 Equations and assumptions 

The reacting Euler, or Zeldovitch-von Neumman-Doering (ZND) equations commonly used 
to model combustion, expressed in Lagrangian coordinates, are 



(2.1) 



d t r - d x u = 0, 
d t u + d x p = 0, 
d t E + d x {pu) = qkc/)(T)z, 
8 t z = -kcf>(T)z, 



2 See problem 29, p. 104 of [CL], cited in [Erl]. 
These estimates, valid on (—00, —Clog |A|] for frequency |A| — > oo, do not seem sufficient for our purposes. 



where r > denotes specific volume, u velocity, E = e + -^u 2 specific gas-dynamical energy, 
e > specific internal energy, and < z < 1 mass fraction of the reactant. Here, k > 
measures reaction rate and q heat release of the reaction, with q > corresponding to 
an exothermic reaction and q < to an endothermic reaction, while T = T(t, e,z) > 
represents temperature and p = p(r, e, z) pressure. 

The equations are of quasilinear hyperbolic type provided that (but not only when) 

(2.2) (p,T) = (p,T)(r,e) and p,p T ,T,T e > 0. 

For simplicity, we assume throughout this paper an ideal gas equation of state and Arrhenius- 
type ignition function, 



s 



(2.3) p = Tr~ l e, T = c -i e, cf>(T) = e" 

where E = e + u 2 /2 is specific (gas-dynamical) energy, c > is the specific heat constant, 
r > is the Gruneisen constant, and £ > is activation energy. Our results on small 
heat-release and high-frequency stability clearly extend to the general case (2.2); however, 
our main results, on the high-overdrive limit, depend in an essential way on the invariances 
associated with (2.3). 

3 Detonation profiles and parametrization 

A right-going strong detonation wave is a traveling-wave solution 
(3.1) 

(u, z)(x,t) = (u, z)(x — st), lim (u, z)(x) = (u-,0), (u, z)(x) = («+, 1) for x > 

of (2.1) with speed s > connecting a burned state on the left to an unburned state on 
the right, smooth for x < 0, with a Lax-type gas-dynamical shock at x = 0, known as the 
Neumann shock. 
Rescaling 



f T + sx T + S 2 t T U T \ 

{x,t,s,T,u,T) -)• I— — , — — ,1, — , ,-2-g J, (z,q,k,£) 

V L L t+ r + s T7_s z / 



'T + sx r + s 2 t 1 r u T \ / g LA; <S 

) J — 2 9 J' \ZiQi K,C ) >■ I Z, 2 o > oi o~ 9 

r + r + s r^s^/ V t+s z t + s z t+s z 

following [Zl], we may take without loss of generality s = 1, r+ = 1 and (by translation 
invariance in u), u + = 0, leaving e + > as the parameter determining the Neumann shock. 
By explicit computation ([Zl], Appendix C), we have then 

(*n - 1 - - r(Te+ + l-f) 

(3.2) u = 1 — r, e - 



r 

(3-3) 

(r + i)(r e+ + 1) - V(r + i) 2 (r e+ + 1) 2 - r(r + 2)(i + 2(r + i) e+ - 2q(z - 1)) 

r = . 

r + 2 



where 

(r + i) 2 (r e+ + 1) 2 - r(r + 2)(i + 2(r + i)e H 



(3.4) < q < q cj :-- 



2T(T + 2) 



and z' = k(f)(c 1 e(z))z; in the simplest case £ = 0, z = e kx . 

The jump at the preceding "Neumann shock" at x = is given (see [Zl], App. C) by 



(3.5) [W] := (l - f (0), f (0) - 1, e+ - e(0), o) 



r 



Taking finally k = 1 by a simultaneous rescaling of x and t if necessary, we can 
parametrize all possible ZND profiles by 

(3.6) (e+,q,£,T), 

where < e + < f7rTTT' < q < q c j(e + ), < £ < oo, and < T < oo. See [Zl] for further 
details. 

3.1 The high- overdrive limit and the scaling of Erpenbeck 

A similar scaling was used by Erpenbeck in [Er3], but with e+ held fixed instead of wave 
speed s. Converting from Erpenbeck's to our scaling amounts to rescaling the wave speed, 
so that T —7- T/s 2 and £ — > £/s 2 , and t — > ts 2 (u is translation invariant, so irrelevant). 
Thus, as noted in [Zl], the high-overdrive limit discussed in [Er3], in which s — > oo with 
u + held fixed, corresponds in our scaling to taking £ = £o e +> Q = 9o e +> an d varying e+ 
from e+ = e c j(qo) (s minimum) to (s = oo), where e c j is determined implicitly by the 
relation q c j(e c j) = qoe c j: that is, the simultaneous zero heat release, zero activation energy, 
and strong shock limit e+ — > 0, £ — > 0, and q — > 0. 

4 Asymptotic ODE theory 

4.1 The conjugation lemma 

Consider a general first-order system W' = A p (x,X)W , W £ C^, A G C, x G R, with 
asymptotic limit A^_ as x — > — oo, where p €. M m denote model parameters and ' d/dx. 

Lemma 4.1 ([MeZl, PZ]). Suppose for fixed 9 > and C > that 

(4.1) \A p - A p _\{x,\) <Ce~ e W 

for x < uniformly for (\,p) in a neighborhood of (Ao,j>o) an d that A varies analytically 
in A and continuously in p as a function into L°°{x). Then, there exists in a neighborhood 
of (Ao,Po) an invertible linear transformation P p {x,\) = I + ^ p {x, A) defined on x < 0, 
analytic in A and continuous in p as a function into L°°[0, ±oo), such that 

(4.2) |^ p | < Cie^>l for x < 0, 



for any < 9 < 9, some C\ = Ci(9,9) > 0, and the change of coordinates W =: P P Z 
reduces W = A P W to the constant- coefficient limiting system Z' = A P _Z for x < 0. 

Proof. See the proof of Lemma 2.5, [Z4], or Lemma A.l, [Zl]. □ 

4.2 The tracking lemma 

Consider an approximately block-diagonal system 

(4.3) W = (^ ^) (x,p)W + 5(x,p)Q(x,p)W, 

where G is a uniformly bounded matrix, S(x) scalar, and p a vector of parameters, satisfying 
a pointwise spectral gap condition 

(4.4) mino-(KMf) - maxcr^Mf) > rj(x) > Ofor all x. 

(Here as usual 'OIN := ^(-^ + N*) denotes the "real", or symmetric part of N.) 

Lemma 4.2 ([MaZ3, PZ, Zl]). Consider a system (4.3) under the gap assumption (4.4), 
with Q p uniformly bounded and rj £ L\ . If sup(5/r))(x) is sufficiently small, then there 
exists a unique linear transformation <3?(x,p) ; possessing the same regularity with respect to 
p as do coefficients Mj and 5@ (as functions into L°°{x)), for which the graph {(Z±, <&Z\)} 
is invariant under (4.3), and 

(4.5) |$ p 0r)| <C f X e f y- v{z)dz S{y)dy < sup {8/rj). 

J —OO ( — CSD,x] 

Proof. See the proof of Lemma A. 4 together with Remark A. 6 in [Zl]. □ 

4.3 A variable-coefficient conjugation lemma 

The key new technical contribution of this paper at the level of asymptotic ODE is the 
following simple observation. Consider a first-order system 

(4.6) W = A p (x;e)W :=M p (x;e)W + e p (x;e)W, x < 0, 
W £ C N , x £ R, p £ R m , with distinguished parameter e — > 0, satisfying 

(4.7) \& p {x,e)\ <Ce 2 e- d£ M, 

(4.8) |»AfP(x,e)| < e5 p (e) + Cee" de \ x \ 
for some uniform C, 9 > 0, all x < 0, where KM := \{M + M*). 



Lemma 4.3. For S p (e) < <5* sufficiently small, and e > sufficiently small, there exists an 
invertible linear transformation P p (x,e) = I + ^ p (x,e) defined on x < such that 

(4.9) \^ p \ < dee- e£lxl/2 for x < 0, 

and the change of coordinates W =: P P Z reduces W = A P W to Z' = M P Z . 

Proof. Equivalently, we construct a solution P p of the (matrix- valued) homological equation 

(4.10) P' = M p P + Q p P, MP:=M p P-PM p , 

satisfying P p — > I as x — > — oo, or, equivalently, a solution *$> p of the integral fixed-point 
equation 

(4.11) T*(x) = f ^ x &(y)(I + V(y))dy, 

J —oo 

where py~^ x is the solution operator of P' = MPP from y to x. 

Denoting by (P : Q) := Trace(P*Q) the Frobenius inner product, and ||P|| := (P, P) 1 ' 2 
the Frobenius matrix norm, we find that 

-||P|| 2 = K(P : P') = 3?(P : M P P - PM P ) = (P : (KM P )P - P(ftM p )) 

<2PM P ||||P|| 2 , 

yielding by (4.8) the bound 

(4.13) 11-^1 < Ce 2eSP{e){x - y) < Ce 2e& " {x - y) . 

For 5* < 9/4 and e > sufficiently small, this implies that T is a contraction on 
L°°(— oo,0]. For, applying (4.7), we have 
(4.14) 

|T*i - r* 2 | (x ) < Ce 2 |^i - ^ 2 |oo f e ee( ~ x - v ^ 2 e eey dy < Cie|*i - * 2 |ooe- 9£|:E|/2 , 

J-oo 

which for e sufficiently small is less than ^l^i — ^loo- 

By iteration, we thus obtain a solution ^ G L°°(— oo,0] of ^ = T^. Further, taking 
^i = \£, ^ = in (4.14), we obtain, using contraction together with the final inequality 
in (4.14), that |^ - T(0)| LO c ( _ oo ^ < ±|* - 0| LO o(_ OOia , ) , yielding, as claimed, \Vl°°(-oo,x)\ < 
2\T(0) Lao{ _ oo>x) \<2C 1 ee- e ^\/2. D 

5 The Evans— Lopatinski determinant 

We now briefly recall the linearized stability theory of [Erl, JLW, Zl, HuZ2]. Shifting 
to coordinates x = x — st moving with the background Neumann shock, write (2.1) as 

6 



W t + F(W) X = R(W), where 



(5.1) 



W 



u 

E 

w 



F 



( —U — ST \ 

Te/r — su 

uTe/r — sE 

—sz 



V 



/ 



\ 



R 







qkz(j)(u) 

\—kz(f>(u) / 



To investigate solutions in the vicinity of a discontinuous detonation profile, we postulate 
existence of a single shock discontinuity at location X(t), and reduce to a fixed-boundary 
problem by the change of variables x — > x— X it). In these coordinates, the problem becomes 
W t + (F(W) - X'(t)W) x = R(W), x/0, with jump condition X'(t)[W] - [F{W)\ = 0, 
[h(x, t)] := /i(0 + , t) — /i(0 _ , t) as usual denoting jump across the discontinuity at x = 0. 

5.1 Linearization/reduction to homogeneous form 

In moving coordinates, W° is a standing detonation, hence (W°,X) = (W°,0) is a steady 
solution of the nonlinear equations. Linearizing about (W°,0), we obtain the linearized 
equations (W t -X'(t)(W Y(x)) + (AW) x = EW, with jump condition X'(t)[W°}-[AW] = 
at x = 0, where A := (d/dW)F, E := {d/dW)R. Computing, we have 

(5.2) 

°\ 





/ o 
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Te 
uTe 



-I 

Vu 





r 



V 



F(e-u 2 ) 
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o 



f 






0/ 
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qkd<t>{T)uz 
c 
kd<f>(T)uz 
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qkd<j>(T)z 

kd<p(f)z 
c 


qk<t>(f) 

-k<P(f)J 



Reversing the original transformation to linear order, following [JLW], by the change of 
variables W —> W — X(t)(W°)'(x), and noting that x-differentiation of the steady profile 
equation F(W°) X = R(W°) gives {A(W°){W )'{x)) x = E(W°)(W°)'(x), we obtain mod- 
ified, homogeneous interior equations Wt + (^4VF)a; = EW together with a modified jump 
condition accounting for front dynamics of X'^fVF ] — L4(VF + X(t)(W°Y)] = 0. 

Lemma 5.1. For < q < q c j, A{x) = dF(W(x)) is invertible for all x, with 
(5.3) \A(x) - A_\ < Ce-"\ x \ 

for all x < 0, some C, rj > 0, where A- := dF(W-). 



defines q c j, marking the boundary of 

□ 



Proof Direct calculation. (The property det A. 
existence of detonation profiles; see [Zl].) 

5.2 The stability determinant 

Seeking normal mode solutions W(x,t) = e xt W(x), X(t) = e xt X, W bounded, of the lin- 
earized homogeneous equations, we are led to the generalized eigenvalue equations (^41^)' = 



(-XI + E)W for x ^ 0, and X(X[W°] - [A(W )'}) - [AW] = 0, where "/" denotes d/dx, or, 
setting Z := AW, to 

(5.4) Z' = GZ, x + 0, 

(5.5) X(X[W°] - [A(W )'}) - [Z] = 0, 
with 

(5.6) G:=(-XI + E)A-\ 

where we are implicitly using the fact that A is invertible, i.e., avoiding the limiting, 
Chapman- Jouget case q = qcj- 

Lemma 5.2 ([Erl, JLW]). For q / qcj, on MA > ; the limiting (n+1) x (n + 1) coefficient 
matrices G± := \im z - ¥ ± 00 G{z) have unstable subspaces of fixed rank: full rank n+1 for G+ 
and rank n for G-. Moreover, these subspaces extend analytically to IRA < — r] < 0. 

Proof. Straightforward calculation using upper-triangular form of G± [Erl, Er2, Zl, JLW]. 

D 

Corollary 5.3 ([Zl, JLW]). For q / q c j, On MA > 0, the only bounded solution of (5.4) 
for x > is the trivial solution W = 0. For x < 0, the bounded solutions consist of an (w)- 
dimensional manifold SpanjZ^ , . . . , Z+}(A, x) of exponentially decaying solutions, analytic 
in A and continuous in parameters (e + ,q,£,T), and tangent as x — > — oo to the subspace of 
exponentially decaying solutions of the limiting, constant- coefficient equations Z' = G-Z; 
moreover, this manifold extends analytically to IRA < — rj < 0. 

Proof. The first observation is immediate, using the fact that G is constant for x > 0, 
with eigenvalues of positive real part. The second follows from standard asymptotic ODE 
theory, Lemma 4.1, Appendix 4, together with the fact that G, by Lemma 5.3, decays 
exponentially to its limit G_ := G(— oo) as x — > — oo, and that G_ by direct calculation 
has a single eigenvalue of negative real part for KA > 0, which extends analytically to 
KA = (by spectral separation from the remaining spectra of G_) and continuously in 
(e + ,q,£,T). □ 

Definition 5.4. We define the Evans-Lopatinski determinant 

L>WA):=det(Zr(A,0), •••, Z~(X,0), X[W°] - [A(W )']) 

= det(Z 1 -(A,0), •••, Z-(X,0), X[W°]+R(W°)(0-)), 

where Z~(X,x) are as in Corollary 5.3. 



The analytic function Dznd is exactly the stability function derived in a different form 
by Erpenbeck [Erl, Er2]. Evidently (by (5.5) combined with Corollary 5.3), A is a general- 
ized eigenvalue/normal mode for RA > if and only if Dznd(X) = 0. 

By duality, the zeros of Dznd agree with those of the ajoint formulation 

(5.8) Dznd(X) = (Z, X[W } + R)\ x=0 , 

where (•, •) denotes complex inner product and Z denotes an analytically chosen solution of 
Z' = —G*Z, x < that is decaying as x — > — oo (see [HuZl, CJLW, Zl]). 

Definition 5.5. For q ^ q c j a ZND detonation is spectrally (or "normal modes") stable if 
the only zero of Dznd (equivalently of Dznd) on 3^A > is a single zero of mulitplicity 
one at A = (necessarily at least multiplicity one by translational invariance). 

6 Continuous dependence and bounded-frequency stability 

Proposition 6.1. For q ^ q c j, Dznd and Dznd ar e analytic in A and vary continuously 
in (e+,q, £,T), uniformly on compact subsets o/{!RA > 0}. 

Proof Immediate, by the construction of the previous section. □ 

Corollary 6.2. Under the parametrization (3.6), ZND detonations are stable with respect 
to bounded frequencies {3ftA > 0} n {|A| < R}, 4 any R > 0, in the small heat release limit 
q < q* sufficiently small for T, 8, e+ bounded and some q* > 0. 

Proof. By continuity with respect to q, and the properties of uniform limits of analytic 
functions, the zeros of Dznd on {3? A > 0} n {|A| < R} converge as q — > to the zeros of 
Dznd with q = 0. Noting that the u and z equations decouple for q = 0, and the profile 
u reduces to a gas-dynamical shock, we find readily from (5.8) that Dznd reduces to the 
Lopatinski determinant for a gas-dynamical shock, which is known (see [Er4, M]) to vanish 
only at a zero of multiplicity one at A = 0. □ 

Corollary 6.3. Under the parametrization (3.6), ZND detonations are stable with respect 
to bounded frequencies {SKA > 0} n {|A| < R}, any R > 0, in Erpenbeck's high-overdrive 
limit, or, in our scaling, q,£,e+ — > 0. 

Proof. Immediate, by Corollary 6.2 and the description of Section 3.1. □ 

7 High-frequency asymptotics and large-frequency stability 

7.1 Approximate diagonalization 

For q bounded away from q c j , k = 1, T bounded, suppressing parameters p = (F,£,q,e + ,X,e), 
rescale x —> x\X\, converting the adjoint eigenvalue ODE Z' = —G*Z to 

Z = B(ex)Z + eC{ex)Z, 



4 That is, Dznd (resp. Dznd) vanishes on this set only at a zero of multiplicity one at A = 0. 



f-l + c 








°\ 





-1-c 














-1 





V o 








-1/ 



where e := |A| _1 , A := A, and B = XA' 1 ' 7 ', C = -(EA~ 1 ) T , with A and £ as in (5.2). 

Noting that A = I J —I, where a is the flux Jacobian for ideal gas dynamics, we find 

from standard gas-dynamical computations that A has distinct eigenvalues — 1, — 1, — 1 ± c, 

where c := -^ — denotes sound speed, and, for < q < qj, these eigenvalues are 

uniformly bounded away from zero as well as from each other. By standard matrix per- 
turbation theory [K], there thus exists a smooth coordinate transformation T = T(B) such 
that T V BT is diagonal. Making the change of coordinates Z = TY, we thus obtain 

(7.1) Y = T~ l (B + eC)TY - T~ X TY =: B X Y + edY, 

where 



(7.2) Bi := T~ X BT = A 



and, since t = (dT/dB)e{dB/dx) = 0{ee- e \ x \) and \E-E(-oo)\ < e~ e \ x \ for x < 0, 9 > 0, 

(7.3) d := T- X CT - e^T^f = (T^CT^-oo) + 0(e- ee ^). 

(* * \ 
, whence 
* hj 

(T~ 1 CT) (— oo) = I T , where E22 = I , ) (dcf)(T)z, 4>(T)) is rank one but (for large 

£ , in particular) not always diagonalizable. Nonetheless, at x = —00, we have 

*,(-») =$)(0,«T_),= ft %%i 

diagonalizable, so, by a further modification of T22, we may take without loss of generality 

(7.4) (T-'CT) 22 (-oo) = (» _ t ; T) ) 

and e(T _1 CT)22 = £ I n uai^ s J +0(ee ' x '). Converting briefly back to x-coordinates 

and applying Lemma 4.1, we thus find that there is a nonsingular coordinate transforma- 
tion W = SX, \S\, IS 1 " 1 ] < C, converting W = eiT^CT^W to the constant-coefficient 

equation X = e f t,A(rp \ ) X. Incorporating this further coordinate change in the 2-2 

block (only), we obtain, finally 

(7.5) X = B 2 X + eC 2 X, 
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(7.6) B 2 = \ 



f-l + c 
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/Q( e -9e\x\ 
y Q ( e -9e\x\ 


0( e -ee\i 

o 2 


V o 





-1- 


-ek<f>(T-)/\) 









where 2 denotes the 2x2 zero matrix. 

Remark 7.1. The reduction just performed, using the conjugation lemma to diagonalize 
the lower righthand block, is a delicate point of the analysis, avoiding potential difficulties 
associated with turning points where E22 becomes nondiagonalizable? 

Using again standard matrix perturbation theory, it follows for e sufficiently small that 
there is a further smooth coordinate transformation 

/ 1 0(ee- fe l £ l) 0(ee- fe l £ l) 0(ee- fe l £ l)\ 

0( ee -<M*|) 1 0( ee -fc|x|j 0(ee- fe l £ l 



Q 



0(ee- 6£ W] 



( £e -fe|£h 

0( ee -fe|*l; 



1 




/ 



such that M := Q 1 {B 2 + eC 2 )Q is block-diagonal, 



(7.7) M 



/A(-l + c) + 0(ee- fe l £ l; 






A(-l 









c) + 0(ee- ee l £ l) 























-A - efc0(T_) ) 



and := — Q l Q = 0(e 2 e ee M) ; taking the equations to the approximate block-diagonal 
form treated in Lemmas 4.2 and 4.3, of 



(7.8) 



w' = mw + ew. 



7.2 High-frequency stability 

Proposition 7.2. Fork = 1, T, £ bounded, and q bounded away from q c j, ZND detonations 
are stable with respect to sufficiently high frequencies; that is, there exists R > independent 
of (T,£,q,e+) such that DzndW 7^ for$t\ > and |A| > R. 

Proof. Case (i) (3JA << e|A|) Equivalently, KA << e, whence, in (7.7), there is a spectral 
gap between the fourth diagonal entry, —A — efc<^(T_), which has real part < — en for 
7] = (f)(T—) > 0, and the first three diagonal entries, A(— 1 ± c) and —A, which have real 
parts > -CKA >> -e. 



'Unimportant for |A| >> 1 + £ , these dominate behavior in the high-activation energy limit £ 
[BZ, Z5]. 
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Applying Lemma 4.2 to (7.8), we find that there is a graph W4 = &(W±, W2, W3) 
<3? = O(e) that is invariant under (7.8), from which we may reduce to an equation 

(7.9) w' = Mw + ew 

on W = (Wi,Ws, W3) alone, with 

'A(-l + c) + 0(ee- fe l £ l) 

(7.1.0) Xr - [j A(-l-c) + 0(ee- e£ l £ l) 

-Ay 

satisfying (4.8) with 5* « e and = 0(e 2 e^ ee ^). 

Applying now Lemma 4.3 to (7.9), we find that there is a coordinate transformation 
W = PX, P = I + 0(e), taking (7.9) to the block-decoupled equation X = MX, of which 
the unique (up to constant multiplier) solution X~ decaying as x — > —00 has value X~(0) 
at x = parallel to (1,0,0, ) T . Untangling coordinate changes, we find that unique (up to 
constant multiplier) solution W~ of (7.8) decaying as x — > —00 is parallel to I + 0(e) times 
(1, 0, 0, 0) T and the unique (up to constant multiplier) solution Z~(0) of adjoint eigenvalue 
equation Z' = —G*Z is parallel to I + 0(e) times the left unstable eigenvector of A _1 (0), 
or (£, 0) T , where £ is the left unstable eigenvector of the gas-dynamical flux Jacobian a(0). 

Likewise, X[W] + R(W(0~)) is parallel to / + 0(e) times [W] = ([u],0) T , whence, 
combining these facts, we find that Dznd(X) = Z~(0) • (X[W] + R(W(0~))) is proportional 
to 1 + O(e) times A (A) := £ ■ X[u], which may be recognized as the Lopatinski determinant 
for an ideal gas-dynamical shock, known by [Er4, M] to be nonvanishing on 3?A > except 
at A = 0, with all constants uniform in model parameters and KA > 0. For |A| sufficiently 
large, therefore, or equivalently, e := |A| _1 sufficiently small, we find that Dznd(X) / 0. 

Case (ii) (5RA > C _1 e|A|) Equivalently, KA > C^ 1 e, whence, applying Lemma 4.2 to 
(7.8), we find that there is a graph (^2,^3, W4) = §(W\), $ = 0(e) that is invariant 
under (7.8), from which we find that the unique unique (up to constant multiplier) solution 
W~ of (7.8) decaying as x — > —00 has value W~(0) at x = parallel to I + 0(e) times 
(1,0,0,0)*. Untangling coordinate changes, and arguing as in the previous case, we thus 
find again that D Z nd(X) = Z~(0) ■ (\[W] + R(W(0~))) is proportional to 1 + 0(e) times 
the gas-dynamical Lopatinski determinant A(A) := £ ■ X[u], £ ■ [u] 7^ 0, hence nonvanishing 
for |A| sufficiently large. □ 

Remark 7.3. Applying Lemmas 4.2 and 4.3 in sequence in this way, one may treat the sit- 
uation arising in the multi-dimensional case (see [Er3] ) of an approximately block-diagonal 
system for which some blocks have a uniform spectral gap and others have uniformly small 
spectral gap. We hope to report on this in future work. 

Remark 7.4. Though we did state it, the arguments above show that there exists a 
change of coordinates Q = l+0(ee ' x ') taking (7.8) to exactly diagonal form W' = MW , 
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where M is as in (7.7). 6 This gives information about the full flow, and not only the 
decaying solution important for the stability theory. Note that, in the exactly diagonal 
coordinates W, the first two entries correspond to coefficients of the first two eigenvectors 
of A in the eigenexpansion of Z, while the second two entries correspond to unknown linear 
combinations of the the coefficients of the third and fourth eigenvectors of A.' 
Likewise, a closer look at the proof reveals the asymptotic description 

(7.11) D ZND (\) = e ClA+c °A(A)(l + O(e)), 

where A(A) := £■ A [it], £■ [it] ^ is the Lopatinski determinant associated with the Neumann 
shock. This can be used as in [HLyZl] as the basis of a convergence study, to obtain practical 
bounds on unstable eigenvalues. Higher order approximants 

Dznd(X) = e ClX+Co+c - lX ~ 1 A(A)(1 + D x e + 0(e 2 )), 

etc., may be obtained by further diagonalizations as detailed in [MaZ3]. 

Corollary 7.5. Under the parametrization (3.6), ZND detonations are stable in both the 
small-heat release and (with Erpenbeck's scaling, Section 3.1) high- overdrive limits. 

Acknowledgement. Thanks to the University of Paris 13 for their hospitality during a 
visit in which this work was partly carried out. Thanks to Olivier Lafitte, Benjamin Texier, 
and Mark Williams for stimulating discussions regarding stability of ZND detonations. 

A Numerical implementation 

For purpose of numerical applications, we provide also a simpler high-frequency argument 
requiring the weaker estimate 



(A.l) $ 



/ 0{ee- e£ \*\) 0(e 2 e- d£2 M) 0(e 2 e- d£ W) 0(e 2 e- e£2 ^)\ 

( e 2 e -fc|x|) o(ee -fe l*l) 0{ee- d£ \ £ \) 0(ee- d£ ^) 

O(e 2 e- 0£ l £ l) 0(ee- e£ W) 0(ee- d£ W) 0(ee- d£ ^) 

\0(e 2 e- e£ ^) 0(ee- e£ W) 0(ee- d£ W) 0(ee- d£ M) J 



rather than O = 0(s 2 e~ £ ' x ') in (7.8), removing the need for the intermediate coordinate 
transformation S in the lower righthand 2x2 block. This avoids an abstract conjugation step 
that is difficult to estimate efficiently numerically, and also provides a slightly simpler proof 
treating all frequencies at once instead of dividing into cases. On the other hand, it provides 
information only about the single mode Z~ associated with the stability determinant, and 
not the entire flow of the adjoint eigenvalue ODE, which may be of interest in more general 
situations. 



6 This follows by separating off scalar diagonal entries with spectral gap from other entries using Lemma 
4.2, to obtain scalar equations w' — mw + 9w with \6\ — 0(e 2 e~ 6,£ ' :c ') for which error 9 can be shown to be 
negligible by explicit exponentiation. (Conjugation of the nonscalar block (7.9) has already been shown.) 

7 Recall, these depend on the abstract conjugation prescribed in going from (7.2) to (7.6). 
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A.l Variable-coefficient gap lemma 

Consider a first-order system (4.6) satisfying (4.7) and 

(A.2) ^M p (x, e) < eS p {e) + Cee- 66 ^ 

for some uniform C,0 > 0, all x < 0, KM := \{M + M*). (Note, in contrast with (4.8), 
that this is a bound from above only.) Assume, further, that there exists a smooth vector 
Vi(x), 1/C < |V?| < C, for which M P V P = 0. 

Lemma A.l. Assuming (A.2), for 5 p (e) < 5* sufficiently small, and e > sufficiently 
small, there exists a solution V p (x,e) of (4.6) defined on x < such that 

(A.3) \{V P - V p ){x)\ < dee- 0e \ x ^ 2 \VV\ for x < 0. 

Proof We seek, equivalently, a solution F p of the integral fixed-point equation 

(A.4) TV(x) = V p (x) + f ^^e(y)F(y)dy, 

J —oo 

where J^v^ x is the solution operator of V = M P V from y to x. From (A.2), we obtain by 
an easy energy estimate like that of (4.12) the bound 

(A.5) 11^1 < Ce eSP{£){x - y) < Ce £S * (x - y) for x > y. 

For 5* < 6/2 and e > sufficiently small, this implies that T is a contraction on 
L°°(— oo,0]. For, applying (4.7) and (A.5), we have 

(A.6) |7Vi - TF 2 | (cc) < Ce 2 !^ - V^ f e e **-M 2 e 9e » dy < C l£ |Vi - F 2 |oce- feN/2 , 

J —oo 

which for e sufficiently small is less than \\V\ — V^loo- By iteration, we thus obtain a 
solution V G L°°(-oo, 0] of V = TV. Further, taking V 1 = V,V 2 = in (4.14), we obtain, 
using contraction together with the final inequality in (A.6), that \V — V^^aof^ x \ = 
\ v ~ T(0)lx,°°(-oo,s) < \\ v ~ 0|i,°°(-oo,s)> yielding (A.3) as claimed. D 

A.2 Alternate high-frequency analysis 

Alternate proof of Prop. 7.2. By inspection, M := (M + 6) - (M + Q)nl satisfies (A.2) 
for M as in (7.7), e := |A| _1 , and 







f@n \ 

o e 22 e 23 e 24 

o e 32 e 33 e 34 

\ o e 42 e 43 e 44 y 



with as in (A.l), with M.V* = for V* := (1,0,0, 0) T , whence, applying Lemma A.l, we 
obtain a decaying solution W(x;e) = ui(x)e x ^~ 1+c ' x V(x;e) of (7.8) converging as 0(e) in 
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e-i-oc ( £e y > o is uniformly bounded above 



relative error to (1, 0, 0, 0) T , where lo(x 
and below. Untangling coordinate changes, and noting that X[W] +R(W(0~)) is parallel to 
l + 0(e) times [W] = ([u],0) T , we thus find that D ZND {\) = Z~(0) ■ (\[W} + R(W(0~))) is 
proportional to 1 + O(e) times the gas-dynamical Lopatinski determinant A(A) := £ ■ X[u], 
£ ■ [u] 7^ 0, hence nonvanishing for |A| sufficiently large. □ 
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